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Background: It is assumed that the introduction of stochastic in mathematical model makes it more adequate. 
But there is virtually no methods of coordinated (depended on structure of the system) stochastic introduction 
into deterministic models. Authors have improved the method of stochastic models construction for the class of 
one-step processes and illustrated by models of population dynamics. Population dynamics was chosen for study 
because its deterministic models were sufficiently well explored that allows to compare the results with already 
known ones. 

Purpose: To optimize the models creation as much as possible some routine operations should be automated. In 
this case, the process of drawing up the model equations can be algorithmized and implemented in the computer 
algebra system. Furthermore, on the basis of these results a set of programs for numerical experiment can be 
obtained . 

Method: The computer algebra system Axiom is used for analytical calculations implementation. To perform the 
numerical experiment FORTRAN and Julia languages are used. The method Runge-Kutta method for stochastic 
differential equations is used as numerical method. 

Results: The program compex for creating stochastic one-step processes models is constructed. Its application 
is illustrated by the predator-prey population dynamic system. 

Conclusions: Computer algebra systems are very convenient for the purposes of rapid prototyping in mathemat¬ 
ical models design and analysis. 

Keywords; stochastic differential equations; “predator-prey” model; master equation; Fokker—Planck equa¬ 
tion; computer algebra software; Axiom system 


I. INTRODUCTION 

This work corresponds our research on mathematical 
models stochastization. This item is interesting due to 
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the following problems: the construction of population 
models from first principles and the introduction of the 
stochastic into such models (the population dynamics is 
studied because of similar models introduction in other 
areas). 

The problem of stochastic term introduction arises dur¬ 
ing mathematical models stochastization. There are sev¬ 
eral ways to solve this problem. The easiest option is an 
in the deterministic equation. But when additive stochas¬ 
tic term is introduced some free parameters that require 
further dehnition appears. Furthermore, these stochas¬ 
tic terms usually interpreted as an external (rather than 
structural) random impact. In this regard, we used 
and improved the stochastic one-step processes models 
construction method, based on master equation UM- 
Stochastic differential equation is considered as its ap- 
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proximate form. It allows to get the model equations 
from general principles. Furthermore, deterministic and 
stochastic parts are derived from the one equation so we 
can regard it as stochastic and deterministic parts con¬ 
sistency. . 

The aim of this work is the software complex devel¬ 
opment for rapid prototyping construction of stochastic 
one-step processes models. This complex consists of two 
blocks. The first block generates the equations of dy¬ 
namic stochastic process model on the principles similar 
to chemical kinetic relations describing the investigated 
process . This block is implemented by means of the com¬ 
puter algebra system - system FriCAS, which is offshoot 
of Axiom. 

The second block is used for the numerical analysis of 
the resulting model. For numerical solution of determin¬ 
istic and stochastic models equations some Runge-Kutta 
different orders methods [^, are used. 

To illustrate the developed system the well-known pop¬ 
ulation model predator-prey is used [m. 

The structure of the paper is as follows. The basic no¬ 
tation and conventions are introduced in Section H Sec¬ 
tion m is devoted to brief introduction to the one-step 
processes stochastization method. Further, in the Sec¬ 
tion m the model under investigation is described. In 
the subsection IIV Al there is a brief reference to standard 
(deterministic) approach, and in the subsection IIV Bl the 
stochastic extension of our model with the help of the 
one-step processes stochastization method is obtained. 

In the Section IV Al we justify selection of the system, 
which implements the model equations generating unit. 
The actual interface of this part of the program complex 
is described in the Section IV Bl . 

The possibility of applying Runge-Kutta methods for 
the analysis of stochastic differential equations is con¬ 
sidered in the Section |vD The software interface of the 
model equations numerical analysis unit is also described 
in this section. Calculations example is based on the 
predator-prey model. 


II. NOTATIONS AND CONVENTIONS 

1. We use abstract indices notation In this nota¬ 
tion tensor as a whole object is denoted just as an 
index (e.g., x*), components are denoted by under¬ 
lined index (e.g., x-). 

2. We will adhere to the following agreements. Latin 
indices of the middle of the alphabet (i, j, k) will ap¬ 
ply to the space of the system state vectors. Latin 
indices from the beginning of the alphabet (a) will 
relate to the Wiener process space. Latin indices 
from the end of the alphabet {p, q) will refer to 
the indices of the Runge-Kutta method. Greek in¬ 
dices (a) will set a number of different interactions 
in kinetic equations. 


3. A Dot over a symbol denotes differentiation with 
respect to time. 

4. The comma in the index denotes partial derivative 
with respect to corresponding coordinate. 


III. ONE-STEP PROCESSES MODELING 

Let’s briefly review the method of one-step processes 
stochastization on the basis of Q. 

We understand one-step processes as Markov processes 
with continuous time with values in the domain of inte¬ 
gers, which transition matrix allows only transitions be¬ 
tween neighbouring portions. Also, these processes are 
known as birth-and-death processes. 

One-step processes are subject to the following condi¬ 
tions: 

1. If at the moment t the system is in state i S 
Z^ 0 )tlien the probability of transition to state i + 1 
in time interval [t, t + At] is equal to At -|- o(At). 

2. If at time moment t the system is in state i S 
Z+,then the probability of transition to state i— \ in 
the time interval [t, t-|-At] is equal to k~At + o{At). 

3. The probability of transition to a state other than 
the neighbouring is equal to o(At). 

4. The probability to remain in the same state is equal 
to 1 — (fc+ -I- k~)At + o(At). 

5. State t = 0 is an absorbing boundary. 

The idea of the one-step processes stochastization 
method is as follows. Based on the patterns of interaction 
we construct a master kinetic equation, expand it into a 
series, leaving only the terms up to and including the 
second derivative. The resulting equation is the Fokker- 
Planck equation. In order to get more convenient model 
we record corresponding Langevin equation. In fact, as 
we shall see , from the patterns of interaction we will 
immediately obtain the coefficients of the Fokker-Planck 
equation (and accordingly, the Langevin equation), so for 
practical use of the method there is no need to construct 
the master kinetic equation. 


A. Interaction schemes 

We will describe the state of the system by a state 
vector cc* inW^, where n is the system dimension (the 
state vector is considered as the set of mathematical 
values, fully describing system). The operator n* S 
X setZ^Q defines the state of the system before the 
interaction and the operator m* G Z"g x Z" g — after.the 
interaction The result of interaction is a system transition 
to another state 0,13. 
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There are s kinds of different interactions that may 
happen in the system, s G So, instead of n® and m® 
let’s consider the operators n®“ G x x Z^q and 
m®“ G Z^Q X Z^Q X Z^Q. 

System elements interaction will be described with 
the interaction schemes similar to chemical kinetic 
schemes 


Q: — l 




+ [s^(x® - r®“,t)p(x® - r®“,t) - s-^{x^)p{x\t)\ } . (5) 

C. Fokker—Planck equation 


ICt 7 V lOL 1 
77, - x-^ ^ m - x% 

k- 


( 1 ) 


Here Greek indexes specify the number of interactions 
and Latin ones specify dimensionality of the system. The 
state change is given by the operator 


r®“ = mf - n®“. 


( 2 ) 


Thus, one step interaction a in forward and opposite 
directions can be written as 

7 7 A 

X —)• X + r-x^ , 

7 7 i 

X —)■ X — r-x -^. 

We can write not in the form of vector equations, 
but in the more traditional form sums: 


UWith the help of the Kramers-Moyal expansion, the 
Fokker-Planck equation la is obtained. For this pur¬ 
pose we will make several assumptions: 

1. there are only small jumps, ie Sa(a:®) is a slowly 
varying function with the change of x®; 

2. p(x®,t) also slowly changes with the change of x®. 
Then in Fokker-Planck equation ® one can shift from 

the point (x® ± r—x ^) to the point x®, and by expanding 
the right-hand side in a Taylor series and dropping terms 
of order higher than the second, we obtain Fokker-Planck 
equation: 

^ = -5, [H®p] + [B^^p] , (6) 

where 


nY'x^Si ^ myx^Si, 

ka 


(3) 


where (5^ = (1,..., 1). 

Also, we will use the following notation: 

n®“ := nf fy, m®“ := (5^ r®“ := rf 5T 


A® := A®(x^^) = r®“ 

H®^ :=S®^■(x^^) = r®“r^“ 


sit - s~ 


st - s" 


a = 1, TO. 


(7) 

As seen from 0, the coefficients of the Fokker-Planck 
equation can be obtained directly from ([5]) and 0, i.e. 
in practical calculations, there is no need to write the 
master equation. 


B. Master equation 


D. Langevin equation 


Transition probabilities per unit of time from the state 
X® to the state x® -I- r-x^ (to the state x® — r-x ^) are 
proportional to the number of x® combination from a set 
of n®“ elements (of x® — combinations from a set of to®“) 
and are given by: 


=«n 


X-! 


^(xi-n®“)!’ 

n 

= kY[ 


xP 

(xl-TOi“)!' 


(4) 


Thus, the general form of the master kinetic equation 
for the states vector x® (it changes by r*—x-^ per step), 
takes the form: 


The Langevin equation corresponds to the Fokker- 
Planck equation: 


dx® = a®df -f 6^dkF“, (8) 

where a® := a®(x^,t), := b\{x^,t), x® G M®® — is the 

system state vector, kF“ G K®®® — TO-dimensional Wiener 
process. Wiener process is implemented as dW = e-\/dt, 
where e ^ Af(0,1) is normal distribution with average 
0 and variance 1. Latin indexes from the middle of the 
alphabet denote the values related to the state vectors 
(dimension of the space is n), and Latin indexes from 
the beginning of the alphabet denote the values related 
to the Wiener process vector (dimension of the space is 
TO ^ n). 

The connection between the equation ([6|) and ([8]) ex¬ 
pressed by the following relationships: 


dp{x\t) 

dt 


A®=a®, 
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We will use Ito interpretation. Under the Ito interpre¬ 
tation, differential of complex functions does not obey 
the standard formulas of analysis. To calculate it rule or 
Ito lemma are used. 

Let / := is a function of a random process 

/ G C^. Then the formula of the differential is [T^ : 


which have the following interpretation. The first rela¬ 
tion means that the prey which eats the food unit im¬ 
mediately reproduced. The second relation describes the 
case when predator absorbs the prey and than it is instan¬ 
taneous reproduced. Only such possibility of prey death 
is considered. Last ratio is a natural predator death. 


d/ = 


5t/ + a7.*+2^a^“/.b 


dt + 6^/.,dlU“, 


where / := a* := a^[x\t), b\ := b\{x\t), and 

dW“ := dfU“(t). 


IV. PREDATOR-PREY MODEL 


A. Deterministic predator—prey model 


All processes are irreversible, so = 0, and 


= k2 


x\ y\ 
{x- 1)! y\ 


X\ 




{x- 1)! {y- 1)! 
x! y\ 


= fciX, 


= hxy, 




Systems with the interaction of two predator-prey pop- 
nlations types are extensively studied and there are a 
lot of varions models for these systems. The very first 
predator-prey model is considered to be a model which 
was obtained independently by A. Lotka and V. Volterra. 
Lotka in [l^ described some hypothetical chemical reac¬ 
tion: 


X % B 

where A, Y are intermediates substances, coefficients 
ki,k 2 ,k 3 are rates of chemical reactions, A is a initial 
reagent, and B is a resultant. As a result was a system 
of differential equations: 


\ X = kix — k 2 xy, 

\y = k2xy - ksy. 

This system is identical to the system of differential 
equations, obtained by Volterra, who considered the 
growth mechanism of two populations with predator- 
prey interaction type. In order to get eqnations @ 
Volterra made a series of idealized assumptions about 
nature of intraspecific and interspecific relationships in 
the predator-prey systenm 


B. Stochastic predator—prey model 

Consider a model of predator-prey system, consisting 
of two individuals species, one of which hunts, second is 
provided with inexhaustible food resources. Let’s intro¬ 
duce the notation, where A is a prey and Y is a predator, 
then we can write the possible processes for the state 
vector xl = (A, Y)^ [3^: 


With the help of the formula ([5]) we have the Fokker- 
Planck equation: 

{A\x,y)p{x,y))+^didj {x,y)p{x,y)) , 

where 

A^ix,y) = s+ix,y) 

B"^{x,y) = s+(x,2/)U“r-^“. 


As a result: 


At(x, y) = kix + ^ k 2 xy + 

0\ j fkix- k 2 xy\ 

-l) '^^y=\k2xy-k3y)' 


2/) = ( Q j (l,0)/cix-k ( ^ \ {-l,l)k2xy + 


-1 


+ (0.-l)^3?/ = 


kix -\- k2xy —k2Xy 
-k2xy k2xy + k^y 


(9) 


In order to write a stochastic differential equation in 
Langevin form ([8]) for predator-prey model, it is enough 
to take the square root of the resulting matrix B*-! in 
Fokker-Planck equation. 



f kxx — k2Xy\ 
\k 2 xy - ksyJ 


dt + b- 




fkix + k2xy 
V -k2xy 


-k2xy \ 
k2Xy + k^yj 


X 

k2 

2A, 

rti = (l,0f 

It should be noted that the specific form of the matrix 
is not written ont because of the extreme awkward¬ 

X + Y 

2Y, 

rl2 = (-l,lf 

ness of the expression. However, with further studies we 

Y 

k3 

0, 

I 

II 

fO 

I 

will need not actually matrix , but its square, i.e. the 
matrix 5^'. 
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V. IMPLEMENTATION OF THE ONE-STEP 

STOCHASTIC PROCESSES MODEL IN THE 
COMPUTER ALGEBRA SYSTEM 

A. Justiflcation of the computer algebra system 
choice 

Let’s consider systems of analytical calculations, Max¬ 
ima and Axiom. Maxima is the first system of analytical 
calculations and it is written in Lisp. Maxima success¬ 
fully runs on all modern operating systems: Windows, 
Linux and UNIX, Mac OS and even on PDA running 
Windows CE/Mobile. Documentation is integrated into 
the program as a handbook with search. There is no dis¬ 
tinction between objects and data in Maxima, and there 
is no clear distinction between the operator and function. 
There is no integrated graphics rendering in the system. 

Unlike Maxima Axiom language is strongly typified 
for better mathematical objects and relationships display. 
The mathematical basis is written in Spad language . Ax¬ 
iom portability is slightly worse: the system runs under 
Linux, UNIX, and graphs does not work under Windows. 
Axiom has its own graphics subsystem. 

In 2007 two Axiom open source forks appeared: Ope- 
nAxiom and FriCAS. Open Axiom is developed by ad¬ 
hering to the ideology of Axiom, problems that occurred 
in the Axiom are eliminated. FriCAS developers reorga¬ 
nized the assembly process, expanded functionality. Fur¬ 
thermore, FriCAS supports not only GCL, which oper¬ 
ates on limited number of platforms, but ECL, Clisp, sbcl 
or openmcl, that allows to run FriCAS under wider range 
of platforms. 

B. Implementation description in the Axiom 
computer algebra system 

Method of one-step processes randomization is orga¬ 
nized as a module for the FriCAS computer algebra sys¬ 
tem. To display all the calculations on the screen the 
variable SH0WCALC;=true is used. To call the method 
you need to use the main function, which has the follow¬ 
ing view: 

osp(Matrix(Integer), Matrix(Integer), 

Vector, Vector, Vector) 

where the hrst argument is before interaction states ma¬ 
trix n*, the second argument is after interaction states 
matrix m®, the third argument is the vector the 
fourth argument is the vector k~, the fifth argument is 
the state vector a;*. Let’s consider the features of the lan¬ 
guage FriCAS on auxiliary functions. For example, the 
function calcProd is used to simplify the calculations 
and s~. In the implementation of the function operator 
of the condition and built-in function reduce are used: 

calcProd : (Matrix(Integer), Vector, Integer, 
Integer) -> Void 


Bufier File Edit Insert Text Paragrai^ Documoit Project Options Help 

osp ([[1.0], [1,1],[0,1]].[[2.0]. [0,2].[0,0]].vector([kl,k2,k3]),vector([0.0, 

0]),v©ctor([x.y])) 

FriCAS will attempt to step throu^ and interpret the code. 

Cannot compile map: calcProd 

We will attempt to interpret the code. 

Compiling function Ai with type (List(Hatrix(lnteger)),Llst( 

Polynomial(Integer)).List(Polynomial(Integer)).NonNegativeInteger 
) -> Matrix(Polynomial(Integer)) 

Compiling function A with type (List(Matrix(Integer)), List( 

Polynomial (Integer)) .List (Polynomial (Integer) ). NonNegativeInteger 
) -> Matrix(Polynomial(Float)) 

Compiling function Bi with type (List(Matrix(Integer)),List( 

Polynomial (Integer)).List(Polynomial(Integer)).NonNegativeInteger 
) -> Matrix(Polynomial(Integer)) 

Compiling function B with type (List(Matrix(Integer)). List( 

Polynomial (Integer)) .List (Polynomial (Integer) ). NonNegativeInteger 
) -> Matrix(Polynaaial(Float)) 

/ —fc2xy + fcla: ^ / k2x.y + klx -k2xy N 
\ —fc3 + fc2xy )\ —fc2xy kiy+k2xy ) 

Figure 1. The output of the module for predator-prey model 
in graphic TgXmacs shell 


calcProd (n, x, a, i) == 

nai:Integer := n(a,i) 

if nai = 0 then 1 else reduce!*, 

[x(i) - j for j in 0..(nai-l)]) 

In the function Bi intermediate calculations for ele¬ 
ments of the matrix are made: 

Bi (rv, sp, sm, i) == rv(i) * 

(trcuispose rv(i)) * (sp(i) + sm(i)) 

In order to nse the module for predator-prey system 
model, we call the function with the following arguments: 

osp ([[1,0],[1,1],[0,1]], [[2,0],[0,2], [0,0]], 
vector([kl,k2,k3]), vector([0,0,0]), 
vector![x,y])) 

Fig. m represents the result obtained in TfjjXmacs shell. 
In fact, we repeated the results obtained in ([S]). 


VI. NUMERICAL EXPERIMENT FOR THE 
PROGRAM COMPLEX 

A. Stochastic Runge—Kutta methods 

Euler-Maruyama method is one of well-known numer¬ 
ical methods for solving SDE, it is a special case of a 
more general Stochastic Runge-Kutta method. Classical 
Runge-Kutta method can be generalized to the case of 
the SDE system (|S]) in the following manner Si: 

XI =xl + hiia\xl ..., AH + ..., xn, 

x\ =xl + hrla\Xl, ...,Xn + • ■ •, 

Indexes k = 1,..., s and I = I,..., n refer to stochas¬ 
tic Runge-Kutta method. J ~ A(0, h) or J ^ y/he, 
e ^ A(0,1) are normal distributed random variables. 
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Such a choice of these numerical values for approxima¬ 
tion is made because the Wiener process is implemented 
as dW = You should also pay attention to double 

summation in the third term of both numerical scheme 
formulas as well as the fact that each number ,..., J” 
should generated separately. 

The method coefficients, as well as for the classical 
analogue, can be grouped into a table called the Butcher 
table: 



Rij 

Rij 


D 

G 


For calculations we used a method with the table 


0 

2/3 

-1 

"T" 


0 0 
0 0 
1 0 
3/4 1/4 


0 

2/3 

-1 

T" 


0 0 
0 0 
1 0 
3/4 1/4 


B. Software implementation description 


Matrix 6 ^ = = '/Wd is calculated numerically 

with the help of the singular value matrix decomposition 
(a subroutine DGESVD from library LAPACK is used). 

For the second subtask scripting language Python was 
chosen (version 3). This language has a wide set of tools 
to work with strings and text files. Except matrices A* 
and additional information about the mathematical 
model was specified as dictionary (standard data type 
in Python), with model name, list of variables, list of 
parameters, initial values of variables, parameters values 
and parameters of the numerical method (integration sec¬ 
tion and step size). 

On the basis of these data, the script automatically 
generates two files functions.f90 and main.f90, where 
the first is a module with functions defining the SDE, and 
the second one is a main program file. While compiling 
these files the third additional module with auxiliary pro¬ 
cedures with Stochastic Runge-Kutta method is added. 


C. The numerical experiment description 

Eor the programs complex work verification a well- 
known predator-prey model was chosen with vector a* 
components 

= ax — f3xy, o? = —yj/ -I- bxy 


The purposes of the programs complex were to auto¬ 
mate the SDE coefficients A* and B'’^ computation with 
the help of general principles described above, and to find 
a numerical solution of the equation obtained by means of 
stochastic Runge-Kutta methods. Erom a programming 
standpoint we can derive three subtasks: 

1 . coefficients A* and B^ generation using the com¬ 
puter algebra system; 

2. generation of source code in languages Eortran and 
Julia, implementing the SDE on the basis of the 
coefficients, saved as a text file; 

3. writing subroutines/fuiictions implementing 
stochastic Runge-Kutta methods in Eortran and 
Julia, and their subsequent compilation together 
with automatically generated source codes. 

As a result of its work Axiom module creates a text file 
which contains the coefficients A* and B'"^ in the following 
form: 

# A 

A[l] 

A[N] 

# B 

B[l,l] B[l,2] . . B[1,N] 

B[N,1] B[N,2] . . B[N,N] 


and matix B'’^: 

ax + jdxy —/dxy 
—fdxy jdxy 'yx ’ 

X is the number of preys, y is the number of predators. 
Coefficients also have the following physical (biological) 
meaning: a is the growth rate of the prey population, 
/3 is a frequency of predators and prey meetings, 7 is 
an intensity of predators death or migration in a lack 
of preys, J is a predator population growth rate on the 
assumption of the excess of the prey. 

During numerical simulations it was taken into account 
that the value of variables a:, y could not be less than zero 
(program stop working when one of the variables becomes 
equal to zero). 

Numerical simulation shows that the addition of 
stochastic to the classical predator-prey model leads to 
the fact that after a certain time death of one of the com¬ 
peting species comes. So, for the following parameters: 
a = 10, /3 = 1, 5, 7 = 8 ,5, J = 1, 8 and the initial values: 
x = 9, 7, x = 6 , 77, victims are first to die, and after 
that predators die due to lack. This case is illustrated 
in Fig. [5] For comparison in Fig. [3] is a graph for the 
deterministic case . 

Under other conditions {a = 10, /3 = 1,5, 7 = 8,5, 
J = 0, 5, X = 22, y = 6, 76) predators die, and the number 
of victims is increasing rapidly, as for their model assumes 
an infinite source of food. Graphics for this case are 
shown in Fig. |T1 and Fig. [3] shown for comparison with 











7 



Figure 2. Stochastic predator-prey model, prey die 



Figure 3. Deterministic predator-prey model 

deterministic case. 


VII. CONCLUSIONS 

This work demonstrates the application of the devel¬ 
oped initial physical system formalization method. The 
system is presented in the form of one or more one-step 
processes. Formalization of the system is done by intro¬ 
ducing the evolution operator. Wherein the analytical 
description of the model requires a lot of routine opera¬ 
tions. To simplify the work we propose to use the com¬ 
puter algebra system (Axiom fork FriCAS). 



Figure 4. Stochastic predator-prey model, predators die 


Figure 5. Deterministic predator-prey model 


We have developed an analytical software package 
block that receives inlet evolution operator and produces 
the SDE, which describes the original model. For numeri¬ 
cal studies of obtained SDE system a second software unit 
that converts the resulting system of equations into the 
program code in fortran and gives its numerical solution 
was developed. Thus, the software system is applicable 
for both analytical and numerical study of the original 
model. 


Currently the software package does not cover all pos¬ 
sibilities, incorporated in the proposed method of for¬ 
malizing the original physical system. Since the original 
system description uses ODE, we should introduce the 
boundary conditions by ties or indicator functions. Par¬ 
tial differential equations can help to solve this problem. 
Further objective is the development of a complete soft¬ 
ware complex for a method of one-step original physical 
system model construction. 
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npeflnocMJiKH: HaniHM KOJijieKTHBOM paspadoTana MeTOflHKa corjiacoBanHoro (saBHcamero ot CTpyKTypti ch- 
CTeMBi) BBefleHHa CTOxacTHKH B fleTepMHHHCTHMecKHe MOflejiH. Ha flaHHOM STane MeTOflHKa orpaHnnena KjiaccoM 
oflHomaroBbix npopeccoB. 

U^ejifa: TIjih onTHMHaapHH padoTbi no coaflannio MOflejieil cjieflyeT aBTOMaTnanpoBaTb Kax mojkho dojinme py- 
THHHbix onepapnii. B pannoM cjiynae npopecc cocTaBjiennH ypaBnennH Mopejin mojkho ajiropHTMHsnpoBaTb n 
peajiHSOBaTb b cncTeMe KOMnniOTepHOH ajire6pbi. KpoMe Toro, na daae sthx pesyjibTaTOB mojkho nojiynnTb n 
nadop nporpaMM pjia npoBepennH MHCjiennoro SKcnepHMenTa. 

MeTopti: 7I,Jia peajinsapnH aHajiHTHnecKHx pacneTOB ncnojibsyeTca cncTeivia KOMObroTepnon ajiredpbi Axiom. 

npoBepenna MHCJienHoro SKcnepnineHTa ncnojibsyiOTca asbixa FORTRAN n Julia. B xanecTBe uncjieHHoro 
MOTopa HcnojibsyeTca MOTop, Pynre-KyTTbi pjia CTOxacTHuecKHx pH4)4)epeHu,HajibHbix ypaBHennii. 

PesyjifaTaTfai: Cospan nporpaMMHbiii KOMnjieKC pjia coapanna CTOxacTnnecKnx MOflejieil opnomaroBbix npopec- 
coB. HpoHJijiMCTpnpoBaHO ero npHinenenne na npnMepe cncTeMbi nonyaapnoHnon pnHaMHKH Tnna «xHm,HHK- 
>KepTBa». ^eTepMHHHCTHuecKHe Mopean paa Taxnx npopeccoB pocTaTonno xopomo nccaepoBanbi, uto nosBoaaoT 
cpaBHHTb noayneHHbie pesyabTaTbi c yjxe nsBecTHbiMn. 

BbiBopfai: CncTeMbi KOMnbiOTepHOH aaredpbi onenb ypodnbi paa peaeit dbiCTporo npoTOTnonpoBanna npn cospa- 
HHH H HCcaepoBaHHH MaTeMaTHuecKHx Mopeaeii. 

Keywords: CTOxacTHaecKHe pH4)4>epeHpHajibHbie ypaBHenHa; Mopeat «xhiii,hhk—?K epTBa»; ocHOBHoe KHnexH- 
aecKoe ypaBHenHa; ypaBHeHHe OoKKepa—riaaHKa; CHCxeMBi KOMHBiOTepHOH aaredpBi; CHCxeMa Axiom 


I. BBEflEHHE 


* eg.eferina@gmail.com 
1 avkorolkova@gmail.com 
t mngevorkyan@sci.pfu.edu.ru 
§ yamadharma@gmail.com 
^ leonid.sevast@gmail.com 

** OnydaHKOBaHO b: One-Step Stochastic Processes Simulation 
Software Package / E. G. Eferina, A. V. Korolkova, M. N. 
Gevorkyan, D. S. Kulyabov, L. A. Sevastyanov // Bulletin of 
Peoples’ Friendship University of Russia. Series «Mathematics. 
Information Sciences. Physics». — 2014. — No. 3. — P. 46— 
59. ; HcxopHBie xeKCTBi: https: //bitbucket. org/yamadharma/ 
art ides-2014-axiom-sdu 


,Z[aHHaH padoTa naxopHTca b pycpe npoBopHMbix naMH 
HCCpepoBaHHH no CTOxacTHsapnH MaTeMaTnnecKnx Mope- 
pen. ,HaHHaH TemaTHKa HHxepecna na-aa cpepyioninx npo- 
6peM: nocTpoenne nonypHpnoHHBix Mopepeli na nepBbix 
npHHpnnoB n BBepenne CTOxacTHKn b Mopepn pannoro 
BHpa (cpepyex aaMexHTB, pto odpainpHCt mbi k nony- 
PHpnoHHOH pHHaMHKe noTOMy, pto naMH nccpepoBapncB 
cxoptne MopepH b ppyrnx odpacTPx). 

Xlpn CTOxacTHaapnn MaxeMaTHPecKnx Mopepeft Boann- 
Kaex npodpoMa BBepenna CTOxacTHPecKoro ppena. Cynte- 
CTByex necKOPBKO cnocodoB cpepaib 3 to. CaMBift npo- 
CTOii Bapnanx — appnxHBHoe podaBpeHne cxoxacxnpecKO- 
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ro Hjiena k fleTepMHHHCTHHecKOMy ypaBHennio. O^naKO 
npH TaKOM BBCfleHHH BOSHHKaiOT CBo6oflHbie napaMCT- 
pbi, Tpe6yioiii;He flajiBHefimero onpe^ejieHHa h3 ^onoji- 
HHTejibHbix coo6pa>KeHHH. KpoMe Toro, flanHbie CTOxa- 
CTHHecKHe HjieHbi oGbihho HHTepnpeTHpyK)TCH kbk Bnein- 
nee (a ne CTpyKTypnoe) cjiyaaftHoe B03;i,eHCTBHe. B cbs- 
3H C 3THM MbI HCnOJIb30BajIH H yCOBCpmeHCTBOBajIH Me- 
TO^ nocTpoeHHH CTOxacTHHecKHx MOflejieft o^HomaroBbix 
npopeccoB na ocHOse ocHOBHoro KHneTHnecKoro ypaBne- 
HHH 0,0. CTOxacTHaecKoe flHcJxJjepenpHajibHoe ypaB- 
Henne paccMaTpHBaeTca kbk ero npHSjiHJKeHnaa (JjopMa. 
3to no3BOJiaeT nojiyaHTb MO^ejibHbie ypaBHenna h3 o6- 
in,Hx npHHpHnoB. KpoMe Toro, ;i,eTepMHHHCTHaecKaa h 
CTOxacTHaecKaa aacTH nojiyaaiOTca h 3 o^noro h Toro ace 
ypaBHenna, hto paccMaTpHBaeTca naMH kbk corjiacoBan- 
HOCTb CTOxacTHaecKoii h ;i,eTepMHHHCTHaecKOH aacTeii. 

IJ,ejibio /i,aHHOfi paGoTbi aBjiaeica pa3pa6oTKa npo- 
rpaMMHoro KOMnjieKca fljia 6biCTporo npoTOTanapoBa- 
HHa nocTpoenaa CTOxacTHaecKHx MOflejieit o^HomaroBbix 
npopeccoB. ^aHHbift KOMnjieKC coctoht h 3 flByx Gjiokob. 
IlepBbiH 6 jiok no cooTHomennaM, anajiornaHbiM cootho- 
mennaM xHMHaecKoft KHncTHKH, onHCbiBaioni;HM nccjie- 
^eMbiii npopecc, renepnpyeT ypaBnenna flHHaMnaecKoii 
CTOxacTHaecKOH MOflejin npopecca. 3 tot 6jiok peajinao- 
Ban na cncTeMe KOMnbiOTepHOH ajire6pbi. B KaaecTBe 
CHCTeMbi fljia peajiHBapHH BbiGpana cncTeMa FriCAS — 
OTBeTBjienne cncTeMbi Axiom. 

BTopoH 6 jiok cjiyacHT ^jia ancjieHHoro anajiHBa noay- 
aennoH MO^ejin. ^Jia ancjieHHoro pemenna ypaBnennii 
fleTepMHHHpoBanHoii n CTOxacTnaecKon MOflejieii ncnoab- 
3yK)Tca MCTO^bi Pynre-KyTTbi pa3Hbix nopa^KOB sa. 

^eMOHCTpai],HH paGoTbi pa3pa6oTaHHoro KOMnjieK- 
ca HcnojibByexcH xopoino HsaecTHaH nonyjiHpHOHHaa mo- 
/i,ejib «xHii],HHK-xcepTBa» ia- 

CipyKTypa CTaTbn cjie^yion^aa. B paB^ejie |TT1 BBe;i,e- 
Hbi ocHOBHbie o6o3HaaeHHa n corjiamenna. B paB^ejie uni 
/i,aeTCH KpaxKoe BBe;i,eHHe b mcto;!, CTOxacxHsapHH o/i,ho- 
maroBBix npoLi,eccoB. ^ajiee, b pa3;i,ejie 113 onncbiBaeTca 
HCCJie;i,yeMaa MO^ejib. Ilpn 3 tom b nojDaBjejie llV Al aaex- 
ca KpaTKaa cnpaBKa no cian^apTnoMy (;i,eTepMHHHCTH- 
aecKOMy) no^xo^, a b no^paB^ejie IIV Bl mbi nojiynaeM 
CTOxacTHaecKoe pacmnpeHne cannon MOflean no MeTO.ii,y 
CTOxacTH3an,HH o^nomaroBbix npopeccoB. 

B pa3fleae IV Al o6ocnoBbiBaeTca Bbi6op cncTeMbi ^aa 
peaanBapHH Gaoxa renepapHn ypaBnennii MO.a,eaH. B pa3- 
fleae IV Bl onncbinaeTca coGctbchho HHiepcjaenc paGoTbi c 
3TOH nacTbio nporpaMMHoro KOMnaeKca. 

B pa3.neae IVII paccMaTpHsaeTca BOBMoacnocTb npnMe- 
nenna MCTOflOB Pynre-KyTTbi flaa anaanaa CTOxacTn- 
necKHx flH4)4)epenn,Haabnbix ypaBnennii n onncbinaeTca 
nporpaMMnbiii nniepcjaeiic Gaoxa nncaennoro anaansa 
MOfleabHbix ypaBnennii. IlpnMep pacneTOB ocnonan na 
MOflean «xnni;Hnx-}xepTBa». 


II. OBOSHAHEHHa H COrJIALUEHHil 

1. B paGoTe ncnoab3yeTca noTapna aGcTpaxTHbix nn- 
;i,excoB 0. B .nanHoii noTapnn TenBop xax peaocT- 
nbiii oGbexT oGoBnanaeTca npocTO nn^excoM (na- 
npnMep, cc*), xoMnonenTbi oGoananaiOTca no^nepx- 
nyTbiM nnflexcoM (nanpnMep, x-). 

2. ByfleM npn;i,ep}xnBaTbca cae.nyioni,nx coraamennii. 
JlaTnncxne nn^excni n3 cepe^nnH aacIsaBnia (i, j, 
k) Gy.nyT OTnocHTbca x npocTpancTBy sexTopoB co- 
CToannii cncTeMbi. JlaTnncxne nn^excbi n3 nana- 
aa aaeJjaBnxa (a) Gy^nyx OTHOcnTbca x npocxpan- 
CTBy BHnepoBCxoro npopecca. JlaTnncxne nn^excbi 
H3 xonpa aa4)aBnTa {p, q) Gy^nyx OTHOcnTbca x nn- 
flexcaM MeTO.ua Pynre-KyTTbi. Fpenecxne nn.zi,ex- 
Cbi (a) Gy.nyT saflanaTb xoannecTBO paanbix B3an- 
MO.a,eiiCTBnii b xnncTnnecxnx ypaBnennax. 

3. Tonxoii na.n cnMBoaoM oGoBnanaeTcn fln4)4)epenn,n- 
poBanne no BpcMcnn. 

4. SanHTOii b nn^nexce oGoBnanaeTcn nacTnan npon3- 
BOflnan no cooTBeTCTByionj,eii xoopflnnaTe. 


III. MOAEJIHPOBAHHE OAHOniAEOBblX 
npoii;EccoB 

xpaTxnii o63op MeTO.a,a CT0xacTn3an,nn o.nno- 
maroBbix npopeccoB na 6a3e paGoTbi 0. 

Ho.!], o/i,nomaroBbiMn npopeccaMn mbi 6y.neM nonnMaTb 
MapxoBCxne npopeccBi c nenpepbiBHbiM BpcMeneM, npn- 
HnMaK)ni;ne Bnanennn b oGaacxn peabix nncea, Maxpnpa 
nepexofla xoTopBix .ii,onycxaeT Toabxo nepexoflbi Mejx.ii,y 
coce.ii,nnMn ynacTxaMn. Taxjxe axn npopeccBi nBBecTHbi 
nofl naBBannHMn npopeccoB po}x.ii,eHnH-rnGean. 

O^HomaroBbie npopeccbi no.ii,nnn5noTCH cae^K)ni,nM 
ycaoBnHM. 

1. Ecan B MOMCHT BpcMcnn t cncTCMa naxoflnTCH b 

cocToannn i S to BepoHTHOCTb nepexo.ua b 

cocToanne z + 1 b nniepBaae BpcMcnn [t, t + Ai] 
pasna k'^At + o{At). 

2. Ecan b momcht BpcMcnn t cncTCMa naxo^nTca b 
cocToannn i S to BeposTHOCTb nepexo.ua b co- 
CToanne z — 1 b nniepBaae speMenn [t, t+At] pasna 
k~At + o{At). 

3. BepoHTHOCTb nepexo.ua b cocToanne, OTannnoe ot 
coceflnnx pasna o{At). 

4. BepoaTHOCTb coxpanenna npejxnero cocToanna 
pasna 1 — (fc+ + k~)At + o{At). 

5. CocToanne z = 0 cctb noraoni,aioni;aa rpannpa. 
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MBTOfla CTOxacTH3an,H o^HomaroBbix npon,eccoB 
cocTOHT B cjie;i,yioiii;eM. Ha ocHOBanHH cxeM BsaHMO- 
flefiCTBHH MbI CTpOHM OCHOBHOe KHHeTHHeCKOe ypaBHe- 
HHe, pacKjiaflbiBaeM ero b pa^, ocTaBjiHH tojibko Hjie- 

Hbl ;i ,0 BTOpOH npOHSBOflHOH BKJHOHHTCJIbHO. HojiyHHB- 
meecH ypaBHenne 6 y;i,eT ypaBHenneM OoKKepa-HjianKa. 

nojiyHeHHH 6ojiee npHBbiHHoro BH^a mo^cjih sanH- 
cbiBaeM cooTBeTCTByromee eMy ypaBHenne Jlan^KeBeHa. 
Ha caMOM flejie, xaK mbi yBH^HM, h3 cxbm B3aHMOfleH- 
CTBHH Mbi cpa3y nojiynacM KOSiJicJiHpHeHTbi ypaBHeHHH 
OoKKepa-HjiaHKa (h, cooTBeTCTBenno, ypaBHeHHH Jlan- 
jKeBena), nosTOMy npn npaxTHHecKOM npHMeneHHH mbto- 
fla CTpoHTb ocHOBHoe KHKCTHHecKoe ypaBHCHHe hbt Heo6- 

XO^HMOCTH. 


A. CxeMfai BsaHMOfleiicTBHa 


CoCTOHHHe CHCTeMbI GyflBM OnHCbIBaTb BCKTOpOM co- 
CToaHHH a;* € R", r;i,e n — pa3MepHOCTb CHCTeMbi (no;; 
BeKTopoM cocToaHHH 6yfleM noHHMaTb mhojkoctbo Maie- 
MaTHaecKHx BejiHHHH, nojiHOCTbK) onHCbiBaioin,Hx chctb- 
My). OnepaTop n* G Z^g x Z^g 3a;i,aeT cocToaniie chctb- 
Mbi BsaHMOflelicTBHH, onepaTop m* G Z"g x Z”g — 
nocjie. B pe3yjibTaTe B3aHMOfleHCTBHH nponcxo^HT nepe- 
xo^ CHCTeMbi B /ipyroe cocxoHHHe 0,[i3- 

B CHCTCMe MOJKeT npOHCXOflHTb S BHflOB pa3JIHHHbIX 

B3aHMO^eHCTBHH, r;i,e s G Z+. HosTOMy bmccto onepa- 
TopoB n* H TO* 6y;i,eM paccMaxpHBaTb onepaxopbi n*“ G 

^^0 ^ ^>0 ^ 


H TO*“ G Z”n X Z"n X Z® 


>0 ^ ^ ^^0- 

B3aHMOfleHCTBHe ajieMenTOB CHCxeMbi Gy^eM onncbi- 
Baxb c noMombio cxeM BsaHMoneftcTBHH, nofloGnbiM cxe- 
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MEM XHMHHeCKOH KHHeTHKH 




kt 


( 1 ) 


3;i,ecb rpenecKHe HH^eKCbi 3a/i;aiOT kojikhcctbo BsanMO- 
fleficTBHH, a jiaTHHCKHe — pa3MepHOCTb CHCTCMbi. H3Me- 
HCHHe cocTOHHHH Gy^BT 3aflaBaTbCH onepaTopoM 


r*“ = TOf - n“. 


( 2 ) 


TaKHM oGpa30M, o^hh mar B3aHMo;i,eHCTBHa a b npa- 
MOM H oGpaTHOM HanpaBjieHHHx mojkho 3anHcaTb b BH^e 


X 


X + r-x -’, 


X — r—x-'. 


Mbi TaKJKe MOJKeM 3anHCbiBaTb (HD He B cJlOpMe BCKTOp- 
Hbix ypaBHeHHii, a b BH;i,e Gojiee Tpa^HpHOHHbix cyMM: 




(3) 


r^e (Jj = (1,..., 1). 

TaKHce Mbi GyflCM Hcnojib30BaTb cjie^ioiii,He oGo3HaHe- 


hhh: 


n*“ := 71*“5^, 


:=rn}^5\ 


B. OcHOBHoe KHKeTHHecKoe ypaBHeHHe 


BepoHTHOCTH nepexo^a b eflHHHpy bpcmchh h3 cocto- 
hhhh X B cocTOHHHe X + r~x-' (b cocTOHHHe X — r—x-') 
nponoppHOHajibHbi cootbctctbchho HHCjiy chocoGob bbi- 
Gopa KOMGHHapHH h3 x^ ho n*“ (h 3 a;* no to*“) h onpefle- 

jiHiOTCH BbipanceHHHMH: 


— 


“'•■in 


='=. n 


(a;i-n“)!’ 
1—1 ^ ' 

■II 


x-\ 


(a;i-TOi“)! 

t—1 ^ ' 


(4) 


TaKHM o6pa30M, oGipnii bh^ ochobhofo KHneTHnecKoro 
ypaBHeHHH ^jih seKTopa coctohhhh x*, H3MeHHioni,erocH 

ioL ■ 

maraMH ajihhbi r—x -', npHHHMaeT bh^: 


dp{x'‘,t) 

m 


+ 


a—1 

s^(x* - r*“,i)p(x* - r*“,i) - S^(x*)p(x*,i)] } . (5) 


C. ypaBHeHHe OoKKepa—IIjiaHKa 


/](ajiee, Hcnojib3yH pa3Jio}KeHHe KpaMepca-Moftajia, 
HOJiynaeM ypaBHenne OoKKepa-HjiaHKa [11]. Ajih stofo 
flejiacTCH HecKOJibKO npe^nojiojKeHHH: 

1. HMeioT MecTO TOJibKO Majibie cxaHKH, T.e. Sq,(x*) hb- 
jiaeTCH 4)yHKH,HeH, Me;i,jieHHO H3MeHHioni,eHCH c h3- 
MeHCHHeM X*; 


2. p{x^,t) TaKJKe Me;i,jieHHO H3MeHaeTCH c H3MeHeHHeM 


Tor;i,a MOJKHO BbinojiHHTb c;i,BHr b H3 TOHKH (x* ± 

r—x^) B TOHKy X* H, pa 3 JiO}KHB npaByio nacTb b pa^ Teit- 
aopa H oxGpocHB aaenbi nopa^Ka Bbime BToporo, noay- 
HHM ypaBHeHHe OoKKepa-HjianKa: 

% = -d^ [B^^p] , 


( 6 ) 
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r^e 


:= A^(x^,t) = r‘^ 
:= B^^(x\t) = r*“r^“ 


si - s: 


si — s. 


a = l.m. 


(7) 

KaK BHflHO H3 0, K034)4)Hn,HeHTiji ypaBHeHHH 
OoKKepa-IIjiaHKa mojkho nojiynHTb cpaay h3 © H 
TO ecTB B npaKTHHecKHx pacHexax sanHCbiBaTb ocHOBHoe 
KHHeTHHecKoe ypaBHenne hot neoGxoflHMOCTH. 


D. VpaBHeHMe JlaHJKeBena 


TaKHx CHCTCM cyipecTByeT 6ojibinoe kojikhoctbo pasHO- 
o6pa3Hbix MO^ejieli. CaMoft nepBoft Mo;i,ejibK) «xhiii;hhk- 
}KepTBa» npHHHTO CHHTaxb Mo;i,ejib, nojiynenHyio nesaBH- 
CHMO Apyr OT flpyra A. JIotkoh h B. BojibTeppoli. JIoTKa 
B CBoeii pa6oTe onHCbiBaji neKOTopyio rnnoxeTHHecKyio 

xHMHHCCKyio peaKpHK) (T5 |: 

A^ X % B 

r^e X^Y — npoMCJKyTOHHbie BeipecTBa, K03cj3(|)HLi,HeHTbi 
kl,k 2 ,k 3 — CKOpOCTH XHMHHeCKHX peaKpHH, A — HCXOfl- 

HbiH peareHT, a _B — npo^yKT peaxpHn. B peByjibxaxe 6 bi- 
jia nojiyHCHa CHCxeMa pHcjicjDepeHpHajibHbix ypaBnenHii 
BHpa: 


ypaBHCHHio OoKKepa-IIjiaHKa cooTBCTCTByeT ypaBHC- 
HHe JlaHSKeBena: 


dx^ = Ydt + bldW^, ( 8 ) 

rpe a* := K '■= a:® G K" — bcktop 

COCTOHHHH CHCTeMbI, G M*” — m-MepHblH BHHepOB- 

CKHH npopecc. BHHepoBCKHH npopecc peajiH3yeTca k&k 
dW = e-\/di, rpe e ^ 7V(0,1) — HopMajibHoe pacnpepejie- 
HHe CO cpepHHM 0 h pncnepcHeli 1. JIaTHHCKHMH HnpeK- 
caMH H3 cepepHHbi aji4)aBHTa o6o3HaHaiOTCH BejiHaHHbi, 
OTHOCHipHecH K BCKTopaM cocTOHHHH (paBMcpHOCTb npo- 
CTpancTBa n), a jiaxHHCKHMH HHpeKcaMH h3 nanajia aji- 
cjjRBHTa oGoBHaaaKITCH BCJIHHHHbl, OTHOCHipHeCH K BCKTO- 
py BHHepoBCKoro npopecca (paBMcpnocTb npocxpancTBa 
TO ^ n). 

IIpH 3TOM CBa3b MCJK/iy ypaBHCHHaMH m H (| 8 ]) Bbipa- 
acacTca cjiepyioipHMH cooTHomcHHaMH: 

A^ = a\ B^^ = blly>‘^. 


Mbi 6 ypeM HcnojibsOBaxb HHTepnpeTapnio Hto. B paM- 
Kax HHTepnpexapHH Hto pHcJjcjjepeHpHaji ot cjiojkhoh 
(JjyHKpHH He HOflHHHaeTca CTanpapTHbiM (JaopMyjiaM ana- 
jiHsa. ^Jia ero BbiancjieHHa Hcnojibsyexca npaeujio hjih 
jieMMa Hmo. 

HycTb / := f{x^,t) — cjaynKpaa ot cjiyaaiiHoro npo- 
pecca / G C^. Torpa cjaopMyaa mAcjaepenpHajia 

6ypeT BbirjiapcTb cjiepyiOHi,HM o6pa30M |l2l |: 


d/ = 






rpe / := Y := a%x\t), b\ := h 

:= dlH“(t). 


IV. MO/];EJIb «XHTTT,HHK->KEPTR A » 

A. .ZI^eTepMHHHCTHHecKaa Mopejifa 

«XHII];HHK-5KepTBa» 

CncTeMbi c BsaHMopeficTBHeM pByx bhpob nonyjiapHii 
THna «xHn];HHK-jKepTBa» ninpoKO HCCJiepoBanbi h paa 


\ X = kix — k2xy, 

\y = k 2 xy - ksy. 

Bxa CHCTeMa noanocTbio coBnapaeT c CHCTeMoii ph(|)- 
4)epeHpHaabHbix ypaBnenHii, noayHeHHofi BoabTeppoii, 
KOTopbiii paccMaxpHBaa MexannsM pocxa HHcaennocTH 
pByx nonyaapnii c BsaHMopeiiCTBHeM Tnna «xHni;HHK- 
}KepTBa». ilaa noayaeHHa ypaBnenHii BoabTeppa peaa- 
ex pap HpeaananpoBaHHbix npepnoaoncennii o xapaxxepe 
BnyXpHBHpOBblX H MeJKBHpOBbIX OTHOmeHHH B CHCXeMe 
«XHHI,HHK-}KepTBa» Q. 


B. CTOxacTHnecKaa Mopeat «XHLpHHK-a<epTBa» 


PaccMOxpHM Mopeab cncxeMbi «xHHi,HHK-}KepTBa», co- 
CToam,yK) h3 oco6eii psyx BnpoB, npnneM opnn h3 hhx 
oxoTHTca, BTopoii — oGecnenen HencnepnaeMbiMn nHHi,e- 
BbiMH pecypcaMH. Baepa oGosnanenna X — jKepxBa, Y — 
xHni;HHK, MOJKHO sanHcaTb BOBMoacHbie npopeccbi (|5|) paa 
BeKTopa cocToanna x- = {X, Y)"^ (l4l - [l7l |: 


X 


2 A, 

rii = (l, 0 f 

X + Y 

k2 

2 y, 

= (- 1 , 1 )^ 

Y 


0 , 

ri3 = (o,-lf, 


KOTopbie HMeiOT caepyiomyK) HHTepnpeTapnio. HepBoe 
cooTHomenne oananaeT, nxo acepTBa, KOTopaa cbepaeT 
epHHHpy nHHi,H, neMepaenno penpopypnpyexca. Bxopoe 
cooTHomenne onncbiBaex noraomenne jKepxBbi xHm,HH- 
KOM H MTHOBCHHOe pCnpOpypHpOBaHHe XHHI,HHKa. 3x0 
epHncTBennaa paccMaxpnBaeMaa bobmojkhoctb rnGean 
jKepxBbi. Hocaepnee cooxHomeHne — 3to ecxecTBennaa 
CMepXb XHHI,HHKa. 

Bee npopeccbi neoGpaxHMbi, noaxoMy s“ = 0, a 


s 


+ 

2 


st{x,y) = k 

{x,y) = k 2 -r- 
[x 

^ 


x\ y\ 
^{x- 1 )!?/! 

x\ y\ 


-l)!(y-l)! 
x\ y\ 

^x\ (y- 1)! 


kix, 

= k2xy, 
ksy- 
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BocnojibSOBaBniHCb cjDopMyjioft ([HI) HMeeM ypaBHenne 
OoKKepa-IIjiaHKa: 

= -di {A\x,y)p{x,y)) + ^didj {x,y)p{x,y)) , 
r^e 


A\x,y) = s+{x,y) r*“, 
B'^ix.y) = s+(a;,y)r*“r-^“. 


TaKHM o6pa30M mbi nojiy^Hjin: 


AKx.y) 




k2xy + 


+ 


BA{x,y) = 




/ kix — k2Xy\ 
\k 2 xy -k^y) ’ 


( 1 , 0)kix + 



(-1, l)k2xy + 


/ 0 \. fkix + k2xy -k2xy \ 

(9) 

JXnsi Toro titoGbi sanncaTb CTOxacTHnecKoe flHcJicijepeH- 
pnajiBHoe ypaBHenne b (jDopMe JlaHiKeBena ® jpnsi mo- 
«xHiii;HHK-}KepTBa», ;i, 0 CTaT 0 HH 0 H 3 BjieHB KBa/i;paT- 
HBIH KOpCHB H 3 nOJiyHeHHOfi MaTpHpBI B ypaBHeHHH 
OoKKepa-IIjiaHKa. 


d 



/kix - k2xy\ 
\k2xy - ksy) 


dt + b- 




fkix -\- k2xy 
-k2xy 


-k2xy \ 
k2xy+ k^yj 


Cjie^yeT saMeTHTb, hto KOHKpeTHbiii bh^ MaTpHpbi 6^ 
He BbiHHcaH H3-3a KpafiHeli rpoMoa^KOCTH BbipaxeHHH. 
BnpoHeM, npH flajiBHeiiniHx HCCjie^OBaHHHx HaM nonaflo- 
6htch He codcTBCHHO MaTpnpa a ee KBa^pai, to ecxB 
Maxpnpa B''^. 


V. PEAJIH3ALi;Ha MO TTE .TTM 
O/THO TTTA rORBTX CTOXACTHMECKHX 
nPOLi;ECCOB b chcteme kombbiotephoh 
AJIPEBP bl 

A. 06oCHOBaHHe BBlSopa CHCTeMBI KOMnBKDTepHoii 
ajireSpBi 

PaCCMOXpHM CHCXeMBI anajIHTHHeCKHX BBIHHCJieHHfl: 

Maxima h Axiom. Maxima anjiHexca caMoii nepBoii ch- 
cxeMoii anajiHXHHecKHx BBiHHCJieHHH, HanHcana ona aa 
a3BiKe Lisp. Maxima ycnemno padoxaex na Bcex co- 
BpeMeHHbix onepapHOHHbix CHCxeMax: Windows, Linux 
H UNIX, Mac OS h ^naace na KITK no;; ynpaBjieHneM 
Windows CE/Mobile. ^HoKyMeHxapHa HHxerpHpoBana b 
nporpaMMy b BH^e cnpaBOHHHKa, ocHan^eHHoro cpe^cxBa- 
MH noHCKa. B Maxima Hex pa3flejieHHH na o6x>eKTBi h 
^aHHBie, Hex aeTKoro paBcpaHHaenna Me}K/i,y onepaxo- 


poM H (|)yHKLi,HeH. B CHCTeMe Hex BCxpocHHoft rpacjDHae- 
CKOii OXpHCOBKH. 

B OTjiHHHe OX Maxima, asBiK Axiom CTporo thhhshpo- 
BaH fpasi jiyamero oxodpaacenHa MaxeMaxHaecKHx o6x.eK- 
xoB H B3aHM0CBa3eH. Bca MaxeMaxHuecKaa dasa Hannca- 
Ha Ha H3biKe Spad. IlepeHOCHMOCTb Axiom ayxb xyace: 
CHCxeMa padoTaex no,zi, Linux, UNIX, a no^ Windows ne 
padoxaeT nocxpoenne rpa(|)HKOB. Axiom odjiaflaex cod- 
CTBCHHOH rpa(|)HHeCKOH HOflCHCXeMOH. 

B 2007 ro.ny y Axiom noaBHjiocb ^na cjDopxa c ox- 
KpbixbiM HCxoflHbiM ko;i,om: OpenAxiom h FriCAS. Open 
Axiom pa3padaxbiBaexcH, npH;i,ep}KHBaacb H;i,eojiorHH 
Axiom, ycxpanaroxca npodjieMbi, Koxopnie BCxpeaajiHCb 
B Axiom. Pa3padoxHHKH npoeKxa FriCAS peopraHH30- 
BauH npopecc cdopKH, pacmnpHjiH cjDyHKpnoHaji. KpoMe 
xoro, FriCAS no^epacHBaei ne tojibko GCL, Koxopbift 
padoxaeT na orpaHHuennoM KOjiHaecxBe njiaxcjDopM, ho h 
ECL, Clisp, sbcl hjih openmcl, axo noBBOJiaeT 3anycKaxb 
FriCAS Ha donee mnpoKOM MHoacecTBe njiaxcijopM. 


B. OnHcaHHe peajiH3au,HH b CHCTeMe 
KOMnBKDTepHoii ajiredpBi Axiom 

Mexofl CTOxacxHBapHH o^HomaroBbix npopeccoB opra- 
HH30BaH B BHfle MO,II,yjIH fljia CHCXeMbI KOMHbHDXepHOH aji- 

redpbi FriCAS. ,r[jiH BbiBO^a Bcex BbiancjieHHii na axpan 
HCHOjibBOBana nepeMCHHaa SHOWCALC := true. ^Zl^jia bbi- 
30Ba Mexofla nyacHO BOcnojibBOBaxbca ochobhoh 4)yHKn,H- 
eii, KOTopaa HMeex cjieflyroipHii bh^: 

osp(Matrix(Integer), Matrix(Integer), 

Vector, Vector, Vector) 

r^e nepBbiH apryMenx — Maxpnpa rij cocxoaHHii ^o Bsan- 
MOfleftcxBHa, BTopoli apryMCHT — Maxpnpa m* nocjie Bsa- 
HMOfleliCTBHH, xpexHH apryMeHT — Bexxop k'^, aexBep- 
xbiii apryMenx — bcktop k~ , naxbiH apryMenx — BCKTOp 
cocTOHHHH cc*. PaccMOxpHM ocodeHHOCTH asbiKa FriCAS 
Ha BCHOMoraxejibHbix 4)yHKn,Hax. HanpnMep, 4)yHKn,Ha 
calcProd HcnojibByexca fljia ynpom;eHHa BbiaHCJieHHii 
H s“. B peajiHBapHH 4)yHKLi,HH Hcnojibsyexca onepaxop 
ycjiOBHa H BCxpoeHHaa 4)yHKn,Ha reduce: 

calcProd : (Matrix(Integer), Vector, Integer, 
Integer) -> Void 
calcProd (n, x, a, i) == 

nai:Integer := n(a,i) 

if nai = 0 then 1 else reduce!*, 

[x(i) - j for j in 0..(nai-l)]) 

B 4)yHKn,HH Bi npoHBBOflHTCa npoMejKyxoaHbie Bbianc- 
jieHHa fljia ajiCMCHTOB Maxpapbi B’’^: 

Bi (rv, sp, sm, i) == rv(i) * 

(trcuispose rv(i)) * (sp(i) + sm(i)) 

Uxodbi BOcnojib30Baxbca MO.ii,yjieM fljia pacaexa MO^e- 
JIH CHCXeMbl «XHLLI,HHK-JKepXBa», BbI3bIBaeM cjDyHKpHK) CO 
cjieflyHDipHMH apryMCHxaMH: 
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Buffer File Eklit Insert Text Paragraf^ Docum«it Project Options Help 

cf a H a X, X ia ^ El j a m~ 

T : \ Z i \ )tt®"<->p!-BCgBR 


osp ([[1.0], [1,1],[0.1]].[[2.0], [0.2],[0,0]].vector([kl,k2,k3]),vector([0.0, 
0]).v©ctor([x.y])) 

FriCAS will attempt to step throu^ and interpret the code. 

Cannot compile map: calcProd 

We will attempt to interpret the code. 

Compiling function Ai with type (List(Hatrix(lnteger)),List( 

Polynomial(Integer)).List(PolynomiaKInteger)).NonNegativeInteger 
) -> Matrix(Polynomial(Integer)) 

Compiling function A with type (List(Matrix(Integer)), List( 

Polynomial (Integer)) .List (Polynomial (Integer) ). NonNegativeInteger 
) -> Hatrix(Polynomial(Float)) 

Compiling function Bi with type (List(Matrix(Integer)),List( 

Polynomial(Integer)).List(Polynomial(Integer)).NonNegativeInteger 
) -> Matrix(Polynomial(Integer)) 

Compiling function B with type (List(Matrix(Integer)). List( 

Polynomial (Integer)) .List (Polyn«aial (Integer) ). NonNegativeInteger 
) -> Matrix(Polynaaial(Float)) 

/ -fc2xy + fcla: ^ / k2x.y + klx -k2xy 
\ —fc3 + fc2xy )\ —fc2xy k3y+k2xy ) 


Ri 


Ri 


BbiT^HCJieHHii HaMH 6biji HcnojibsOBan mbto;!, c Ta6- 

jiHuefi 


0 

2/3 

-1 

"T" 


0 0 
0 0 
1 0 
3/4 1/4 


0 

2/3 

-1 

"T" 


0 0 
0 0 
1 0 

3/4 1/4 


Phc. 1. PesyjibTaT pa6oTbi MOflyjia fljia MOflejiH «xHmHHK- 

*epTBa» B rpa4>HMecKoii o6ojiOHKe TlgXmacs ^ Onucauue nporpaMMHofl peajiHsapHH 


osp ([[1,0],[1,1],[0,1]], [[2,0], [0,2], [0,0]], 
vector([kl,k2,k3]), vector( [0,0,0]), 
vector( [x,y])) 

Ha pHC. [T] npe^cxaBjieH peayjibxax, nojiynaromHiicH b 
oGojio’XKe Tj^Xmacs. OaKXHnecKH, mm noBxopnjiH peayjib- 
xaxbi, nojiyTxeHHbie b 


VI. HHCJIEHHblH 3KCnEPHMEHT B 

PAMKAX nPOEPAMMHOrO KOMHJIEKCA 

A. CTOxacTH'iecKHe MeTOflfai Pyare—KyxTfai 

O^HHM H3 xopomo HSBecxHbix HHCJieHHbix MCxoflOB pc- 
meHHH C/I,y HBjiaexcH Mexo^n, Biijiepa-MapyaMbi, koxo- 
pbiH, HBjiaexcH aacxHbiM cjiyaacM 6ojiee o6in,Hx cxoxa- 
cxHHecKHx Mexo;i,OB Pynre-Kyxxbi. KjiaccHHecKHii Mexo;i, 
Pynre-Kyxxbi mohccx 6bixb o6o6meH na cjiynali cncxeMbi 
C/I,y ([8]) cjie^yiomHM oGpaaoM ii: 

XI =xi + hRia^xi, ...,xn + i?x , xn, 

=4 + hria^Xl, ...,Xn+ fLrbliXl, ...,Xn 

Hn^eKCbi k = l,...,svLl = l,...,n oxHOCHxca k cxo- 
xacxHHecKOMy Mexo^y Pynre-Kyxxbi. J ^ -/V(0, h) hjih 
J ~ Vhe, e ^ iV(0,1) — HopMajibHO pacnpe^ejieHHbie 
cjiynaHHbie BejiHHHHbi. Bbi6op xaxnx BejiHHHH pjisi hhc- 
jieHHOH annpoKCHMapHH o6x.HCHHexcH xcm, hxo BHnepoB- 
CKHe npopecc peajiHsyexca kbk dW = ev/di- TaxjKe cjie- 
^ex oGpaxHXb BKHManne na flBoliHoe cyMMHpoBaHne b 
xpexbBM cjiaraeMOM oGchx 4)opMyji ’xhcjichhoh cxeMbi, a 
xaKJKe na xo, hxo Kaxc^oe hhcjio , J" ;],oji}kho re- 

nepHpoBaxbCH ox^ejibHO. 

Ko3(J)4)Hii,HeHXbi Mexo;i,a, xaKJxe xax h ^jih KjiaccHne- 
CKoro anajiora, mojkho crpynnnpoBaxb b xadjinpy, nasbi- 
BaeMyK) ma6AU'u,eu Bamuepa: 


Hpn HanncaHHH KOMnjieKca nporpaMM cxaBHjiacb sa- 
flana xax aBxoMaxHsapHH BbiHHCjienHH KoacjxjDHpHeHxoB 
A* H B*-' C/jy c noMombK) onncbiBaeMbix Bbiine o6- 
Ln,Hx npHHpHnoB, xaK h HHCjienHoe peineHne nojiy’xenHoro 
ypaBHCHKH c noMombK) cxoxacxHHecKHx Mexo^OB Pynre- 
Kyxxbi. C xo^iKH speHHH nporpaMMHpoBanHH sa^ana pac- 
najiacb na xpn no^aa^aHH: 

1. renepapHH KOScJxjDHpHeHxoB A* h npn noMOipH 
CHCxeMbi KOMHbioxepHOH ajireGpbi; 

2. reHepapHH hcxo^hbix xcrcxob nporpaMM na hsm- 
Kax Fortran n Julia, peajinsyromnx C/I,y na ocho- 
Be nojiyneHHbix KOScJxJrHpneHxoB, coxpaneHHbix b 
BH^e xeKCxoBoro dpsAjia; 

3. Hanncanne noflnporpaMM/cjayHKpHii, peajinayio- 
ni;Hx cxoxacxnnecKne Mexo^bi Pynre-Kyxxbi na 
HSbiKax Fortran n Julia, n flajibnenmaa: KOMunjin- 
pHH Hx BMecxe c aBxoMaxnnecKH crenepnpoBannbi- 
MH HCXO/I,HbIMH XeKCXaMH. 

Mo;i,yjib Axiom b pesyjibxaxe CBoeii pa6oxbi co3;i,aex 
xeKCxoBOH cjaanji, b KOxopoM co;i,ep}KaxcH KOScjDcJrHpHen- 
xbi A® H B cjie;i,yioni;eM Bn^e: 

# A 

A[l] 

A[N] 

# B 

B[l,l] B[l,2] . . B[1,N] 

B[N,1] B[N,2] . . B[N,N] 

Maxpnpa 6^ = = 'JB'-i BbinncjinexcH nnc- 

jienno npn noMorpn CHnryjiHpnoro pasjioncennH Maxpn- 
pbi (ncnojibBOBana no/i,nporpaMMa DGESVD GnGjinoxe- 
KH LAPACK). 

pemeHHH Bxopoli noflaa^ann 6biji BbiGpan nnxep- 
npexnpyeMbiH hbmk Python (nepcnn 3). /Hanubin hbmk 
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oGjia^aeT niHpoKHM naGopoM cpeflCTB ^jih paGoTbi co 
CTpOKaMH H TeKCTOBbIMH (JjaftjiaMH. KpOMe MaTpHU, A’^ H 
Gbijia aa^ana ;i,onojiHHTejibHaH HHcJjopMapHH o Maie- 
MaTHa:ecKOH Mo;i,ejiH b BH;i,e cjiOBapa (cTanflapTHbiii thh 
^ aHHbix B Python), r,ne Gbijin yKasaHHbi: nasBaHHe mo- 
^nejiH, cnncoK nepcMeHHbix, chhcok napaMOTpoB, nanajib- 
Hbie SHaneHHH nepeMennbix, SHanenHa napaMOTpoB h na- 
paMCTpbi HHCJieHHoro MOTO^a (oTpeaoK HHTerpnpoBaHHa 
H BCJiHaHHa niara). 

Ha ocHOBe 3 thx flannbix cpenapHii aBTOMaTnaecKH re- 
nepnpyeT ,nBa cjaanjia functions.f90 n main.f90, r,ne 
nepBbiii — 3 to Mo;i,yjib, co,nepjKaiii,HH (jDyHKpHH, 3a,naio- 
ipne C/jy, a BTopoft — ochobhoh (|)aHji nporpaMMbi. Hpn 
KOMHHJiapHH K 3THM (|)aHJiaM floGaBJIHIOTCH eipO Tpn MO- 
^jia CO BcnoMoraTCjibHbiMH npoLi,e,ii,ypaMH, b hhcjic koto- 
pbix co,zi,epjKHTCH CTOxacTHuecKHH MCTO^a, PyHre-KyTTbi. 


C. OnHcaHHe HHCJieHHoro 3KcnepHMeHTa 

B KaaecTBC MO,zi,ejiH fljia BepncjiHKapHH paGoTbi kom- 
njiCKca nporpaMM Gbijia BbiGpana xopoino nsBecTnaa mo- 
flcjib «xHiii;HHK-}KepTBa», 3a,ii,aBaeMaa BCKTopoM a* c kom- 
nOHCHTaMH 

= ax — (3xy, o? = —7?/ + bxy 
H MaTpnpeH : 

ax + ^xy —fixy 
—(5xy fixy + 7a; ’ 

X — HHCJIO JKCpTB, y — HHCJIO XHipHHKOB. KoOcjacJjHpHeH- 
Tbi ace HMeiOT cjie^nyiomHH (JanaHHecKHii (GHOjiornaecKHH) 
CMbicji: a — HHTencHBHOCTb pocTa nonyjiapHH acepTB, 
P — aacTOTa BCTpean xhiii;hhkob h acepTB, 7 — HHTen- 
CHBHOCTb CMepTH HJIH MHrpapHH XHipHHKOB B yCJIOBHH 

He^HOCTaTKa acepTB, 5 — HHTeHCHBHOCTb pocaa nonyjia- 

pHH XHin,HHKOB npH yCJIOBHH HaGbITKa JKOpTB. 

Hpn HHCJieHHOM MOflejinpoBaHHH yuHTbiBajiocb, hto 
anaaeHHe nepeMenHbix x, y ne Moacei Gbitb Menbine nyjia, 
H npH oGpaipeHHH b hojib o^hoh h3 nepeMenHbix paGoTa 
nporpaMMbi npeKpaipajiacb. 

HncjieHHoe MOflejinpoBanne noKaaajio, hto .noGaBjie- 
HHe CTOxacTHKH B KjiaccHHecKyio MO,nejib «xhiii;hhk- 
}KepTBa» npHBOflHT K TOMy, HTO HO nponiecTBHH onpe- 
flejieHHoro bpcmchh HacTynaex CMepxb oflHoro h 3 Konny- 
pHpyioni,Hx bh,hob. Thk, npn cjie;i,yioiii;Hx ananeHHax na- 
paMexpoB: a = 10, /3 = 1, 5, 7 = 8, 5, J = 1,8 n cjieflyio- 
ru,Hx nanajibHbix ananennax: x = 9, 7, x = 6, 77 nepBbiMH 
HornGaiOT acepTBbi, a aa hhmh h xhiii;hhkh bbh^ He.no- 
CTUTKa nnipH. .HaHHbiii cjiynaii npoHjunocTpnpoBaH na 
pHC. m /Ijih cpaBneHHH na pnc. [3] npencTaBjien rpaiJiHK 
nJia nexepMHHHpoBaHHoro cjiynaa. 

Hpn npypHx ycjiOBHHx {a = 10, ,8 = 1,5, 7 = 8,5, 
6 = 0,5, X = 22, y = 6,76) nornGaiox xhiii;hhkh, a 
HHCjieHHOCxb acepxB cxpeMHxejibno Boapacxaex, xax kuk 
njia HHx Monejib npennojiaraex GecKoneanbiH hcxohhhk 



Phc. 2. CxoxacTHMecKaa MOflejib «xHmHHK->KepTBa» — 
CMepxb acepTB 



Phc. 3. H^exepMHHHpoBaHHaa MonejiB «xHni,HHK-acepTBa» 

HHipH. PpaiJiHKH fljin na-HHoro cjiyaaa HaoGpaacenbi na 
pHC. |4j a na pnc. On-na cpaBnenna noKaaan nexepMHHH- 
poBaHHbiii cjiynaii. 


VII. 3AKJIIOHEHHE 

B paGoxe nponeMoncxpHpoBaHO npaMenenne paapaGo- 
xanHoro Mexona 4)opMaaH3an,HH ncxonnon (|)H3HnecKOH 
CHCxeMbi, npencxaBHMOii b BH^e OflHoro hjih necKoabKHx 
onHoniaroBbix npopeccoB. OopMajinaapna chcxcmbi npo- 



Phc. 4. CTOxacTHuecKaa MOjiejib «xHiii,HHK-jKepTBa» — 
CMepXb XHHIHHKOB 






















Phc. 5. TJeTepMHHHpoBaHHaa MOflejib «xHiii,HHK-jKepTBa» 


H3BOflHTCH nyTCM BBe^eHHH onepaTopa 3 bojiioli,hh. IIpH 
3TOM fljiH nojiyneHHH anajinTHHecKoro onncanHa MO^ejiH 
Tpe6yeTCH npHMeHenne 6ojiBinoro nncjia pyTHHHbix one- 
pau,HH. ynpomeHHH pa6oTbi npefljioJKeno ncnojibso- 
BaTb cpeflCTBa KOMUBiOTepHOH ajireGpBi — kjioh Axiom 
FriCAS. 


HaMH 6biji paapaGoTan anajinTHHecKHii nporpaMMHbiii 
6 jiok naxeTa, nojiyaaiomHH na Exo^e onepaTop sbojiio- 
pHH, a na Bbixofle BbiflaiomHH CHCTOMy C^Y, onHCWBaio- 
myio HCxoflHyio Mo;],ejiB. ^jih HHCjienHoro HCCJiCflOBanHa 
nojiyHeHHOii TaxHM o6pa30M cncTeMBi C/JY pa3pa6oTaH 
BTopoii nporpaMMHbiii 6 jiok, KOTopBiii npeo6pa3yeT nojiy- 
Hennyio b FriCAS CHCTOMy ypaBnenHit b nporpaMMHbiii 
KOfl Ha H3biKe fortran h na Bbixo^e ^aex ee nncjiennoe pe- 
menne. TaxnM o6pa30M, pa3pa6oTaHHbiii nporpaMMHbiii 
KOMHjieKC noaBOjiHOT nponecTn oflnoBpeMeHHO anajiHTH- 
HecKoe H HHCJienHoe HCCjie;],OBaHHe hcxo^hoh MOflCjin. 

Ha ;i,aHHbiii momcht paapaGoTannbili nporpaMMHbiii 
KOMHjieKC oxBaTbiBaeT ne Bce bohmohchocth, aajioJKen- 
Hbie B npe;i,jio>KeHHOM MeTo;i,e 4)opMajiH3aH,HH hcxo^hoh 
4)H3HHecK0ii cncTeMbi. Tax b cjiynae onncanHs hcxo^hoh 
CHCTeMbi c HOMOipbio O/JY, rpaHHHHbie ycjiOBna npnxo- 
flHTCH BBOflHTb C nOMOin,bIO CBH3eH HJIH HHflHXaXOpHblX 

(JiyHxpHH. B xojxe BpeMa npaMeHenne flH4)4>epeHH,Hajib- 
Hbix ypaBHCHHii b nacxnbix nponaBO^Hbix Mojxex no3- 
BOjiHXb pemHXb 3xy npo6jieMy. /^ajibneiimHe aa^ann — 
pa3pa6oxxa nojinoro nporpaMMHoro xoMnjiexca npHMe- 
Henna MCxoflHXH nocxpoenna o^nomaroBbix MO^ejinii nc- 
xo^Hoii 4)H3HHecxoii cncxeMbi. 
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